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The superconducting gap structure of heavy fermion UPd2Al3, in which unconventional supercon- 
ductivity coexists with antiferromagnetic (AF) order with atomic size local moments, was investi- 
gated by the thermal conductivity measurements in a magnetic field H rotating in various directions 
relative to the crystal axes. The thermal conductivity displays distinct two-fold oscillation when H 
is rotated in the plane orthogonal to the basal a&-plane, while no oscillation was observed when H 
is rotated within the basal plane. These results provide strong evidence that the gap function A(fe) 
has a single line node orthogonal to the c-axis located at the AF Brillouin zone boundary, while 
A(fc) is isotropic within the basal plane. This gap structure indicates that the pairing interaction in 
neighboring planes strongly dominates over the interaction in the same plane. The determined nodal 
structure is compatible with the resonance peak in the dynamical susceptibility observed in neu- 
tron inelastic scattering experiments. Based on these results, we conclude that the superconducting 
pairing function of UPd2Al3 is most likely to be d-wave with a form A(fe) = Ao cos{kzc). 

PACS numbers: 74.20.Rp, 74.25. Fy, 74.25.Jb, 74.70. Tx 
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I. INTRODUCTION 



The relationship between superconductivity and mag- 
netism is a fundamental problem in the physics of 
strongly correlated electron systems. Over the past two 
decades, unconventional superconductivity with symme- 
tries other than s-wave symmetry has been found in 
several classes of materials, including heavy fermion 
(HF), high- Tc, ruthenate, and organic superconductors. 
Unconventional superconductivity is characterized by 
anisotropic superconducting gap functions belonging to a 
nontrivial representation of the crystal symmetry group, 
which may have zeros (nodes) along certain directions in 
the Brillouin zone The nodal structure is closely 

related to the pairing interaction of the electrons and it 
is widely believed that the presence of nodes in almost all 
unconventional superconductors discovered up till now is 
a signature of a magnetically mediated interaction, in- 
stead of the conventional electron-phonon mediated in- 
teraction. HF compounds containing Ce and U atoms 
are known to exhibit a rich variety of unconventional su- 
perconductivity associated with their peculiar magnetic 
properties In these materials, as the temperature is 
lowered, /-electrons, which are well localized with well- 
defined magnetic moments at high temperatures, begin 
to become delocalized due to the hybridization of atomic 
and conduction electron wave functions. Eventually, at 
very low temperatures, the /-electrons appear to be- 
come itinerant with heavy effective electron mass about 
one hundred times the free electron mass. In HF su- 
perconductors, the strong Coulomb repulsion within the 



atomic /-shells leads to a notable many-body effect and 
often gives rise to Cooper pairing states with angular mo- 
mentum greater than zero. Therefore understanding the 
nodal structure of HF superconductors is important in 
understanding the physics of unconventional supercon- 
ductivity in strongly correlated system. 

Among HF superconductors, UPd2Al3 has aroused 
great interest because it displays quite unique proper- 
ties. In UPd2Al3, superconductivity with heavy mass 
(Sommerfelt specific coefficient 7=140 mJ/K^mol) oc- 
curs at Tc=2.0 K after antiferromagnetic (AF) order- 
ing with atomic size local moments (/i —0.85fiB) occurs 
at TAr=14.3 K Below Tc, superconductivity coex- 
ists with magnetic ordering. The ordered moments are 
coupled ferromagnetically and lie in the basal hexago- 
nal a6-plane along the a-axis. The ferromagnetic sheets 
are stacked antiferromagnetically along the c-axis with 
wave vector Qo=(0,0,7r/c), where c is the c-axis lattice 
constant ^ (For the structure of the hexagonal basal 
plane, see the inset of Fig. 1). Due to a pronounced 
magnetic anisotropy, the aligned spins tend to fluctuate 
in the basal plane. The presence of large local moments is 
in contrast to other HF superconductors, in which static 
magnetic moments are absent or very small if present 
|3,la|- In UPd2Al3, both the superconductivity and AF 
ordering are associated with the 5/ electrons in the U 
atoms. The presence of the two subsystems, partly local- 
ized and partly itinerant, has been extensively studied 
in the context of the dual nature of strongly correlated 
electrons Hllll^]. 

In the superconducting state of UPd2Al3, two notice- 
able features have been reported. The first is the modu- 
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lation of the tunneling conductivity above the supercon- 
ducting gap at an energy of about 1.2 meV in cross- type 
tunnel junctions UPd2Al3-A10a;-Pb, on a UPd2Al3 thin 
film 10]. Furthermore, neutron inelastic scattering ex- 
periments have revealed a strongly damped AF spin- wave 
excitation (magnetic exciton) with an excitation energy 
of about 1.5 meV at and in the vicinity of Qo well below 
T/v. Such an anomaly in the tunneling spectrum is well 
known from ordinary strong coupling electron-phonon su- 
perconductors like Pb. In analogy, the "strong coupling 
anomalies" observed in UPd2Al3 were attributed to the 
strong interaction between the HF quasiparticles and AF 
spin-wave excitations. The second feature is the appear- 
ance of a "resonance peak" in the neutron inelastic scat- 
tering in the vicinity of Qo at about 0.3 meV well below 
Tc U llll ll^ . This peak has been interpreted as the 
creation of quasielectron/quasihole pairs from supercon- 
ducting condensed pairs by magnetic neutron scattering, 
which in turn appears to be closely related to the sign 
of the superconducting order parameter under the tran- 
sition k + Qo It should be noted that a res- 
onance peak in the superconductin g sta te has also been 
reported in d-wave high- Tc cuprates [lj| and borocarbide 
superconductors but its origin is controversial. 

A major outstanding question about UPd2Al3 is the 
nature of the microscopic pairing interaction responsi- 
ble for the superconductivity. To elucidate the pairing 
mechanism, identification of the symmetry of the super- 
conducting order parameter is of primary importance. 
In particular, the gap structure of UPd2Al3 is expected 
to provide most valuable information on the relation- 
ship between the superconductivity and magnetism. The 
gap function of UPd2Al3 has been the subject of ex- 
tensive studies. Thermodynamic and nuclear magnetic 
resonance (NMR) measurements have revealed that the 
gap function is anisotropic JJ, 17]. In particular, the 
NMR relaxation rate Tj~^, which exhibits a behavior 
in the superconducting state over four orders of mag- 
nitude down to O.lTc, along with the absence of the 
Hebel-Slichter coherence peak, indicate the existence of 
line nodes in the gap function |17|. Pauli limiting Hc2 
pl |. NMR Knight shift reduction [l3|, and /iSR reduc- 
tion in spin susceptibility '19] below Tc favors a spin sin- 
glet state. Thus spin-singlet superconducting symmetry 
with line nodes in UPd2Al3 seems to be confirmed. 

It has been pointed out that the resonance peak around 
Qo in the neutron inelastic scattering spectrum implies 
an order parameter displaying sign inversion on transla- 
tion by Qo on the Fermi surface, due to the coherence 
factor in the dynamical susceptibility of the conduction 
electrons. In the simplest scenario, A(fc) = — A(fc + Qo) 
is inferred, which suggests the existence of nodes orthog- 
onal to the c-axis (horizontal node), (12,, 13, 20J. Here 
A(fc) is the superconducting order parameter. Although 
the coherence factor argument may place a constraint 
on the sign of the order parameter, the neutron scat- 
tering intensity mainly arises from a small minor band, 
which seems to be irrelevant in the occurrence of super- 



conductivity, as suggested in Ref. Moreover, the 

neutron scattering results do not provide information on 
the position and number of the horizontal nodes and of 
the gap structure within the basal plane (vertical node). 
The existence of horizontal nodes was also suggested by 
the tunneling conductance along the c-axis 10]. However 
the tunnelling conductance in anisotropic superconduc- 
tors strongly depends on the tunneling process of the 
junction, whether the process is coherent or incoherent, 
and on the transfer matrix originating from the atomic 
orbital at the surface, as extensively studied in high- Tc 
cuprates 0, [S^, HI]. In fact, the typical d-wave gap 
structure has never been observed in c-axis tunneling ex- 
periments of high- Tc cuprates. On the other hand, the 
presence of line nodes perpendicular to the basal plane is 
suggested by the anisotropy of the thermal conductivity 
|24l |. Thus, despite great experimental effort, the detailed 
structure of the superconducting order parameter is still 
an unresolved issue, therefore, there is a need for a probe 
of the quasiparticle excitation with spatial resolution in 
fe-space. 

The purpose of this work is to determine the super- 
conducting order parameter A(fc) of UPd2Al3 by mea- 
suring the angle resolved magneto-thermal transport, 
which has proven to be a powerful probe for determining 
low energy quasiparticle excitations, including direction 

The thermal conductivity is a unique transport quantity 
that does not vanish in the superconducting state, re- 
sponding to unpaired quasiparticles that carry the heat. 
We provide direct evidence that the gap function has 
a horizontal node located at the AF zone boundary 
kz = ±7r/2c, and is isotropic within the basal plane. 
These results place strong constraints on models that at- 
tempt to explain the mechanism of the unconventional 
superconductivity of UPd2Al3. We discuss the pairing 
mechanism infered by the gap structure. 



II. EXPERIMENTAL 

High quality single crystals of UPd2Al3 with Tc — 
2.0 K were grown by the Czochralski pulling method in 
a tetra-arc furnace. The residual resistivity ratio was 
55 along the 6-axis and 40 along the c-axis, indicating 
the crystals were of the highest crystal quality currently 
achievable. The upper critical fields in H parallel to the 
a-, b~, and c-axes at 0.4 K, determined by the resistive 
transitions, were i/c°2=3-2 T ,77c^2=3-l T, and H^2^3.8 T, 
respectively. We measured thermal conductivity along 
the c-axis of the hexagonal crystal structure, k^z (heat 
current q || c) and along the 6-axis Kyy {q |j b) by the 
steady-state method. The sample dimensions for the Kzz 
and Kyy measurements were 0.32x0.34x2.50 mm"^ and 
0.22x0.50x1.80 mm^, respectively. To apply H with 
high accuracy relative to the crystal axes, we used a sys- 
tem with two superconducting magnets generating mag- 
netic fields in two mutually orthogonal directions and a 
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FIG. 1: Temperature dependence of the thermal conductivity 
in zero field in the b-axis Kyy {q\\ b) and c axis n^z {q\\ c). Inset 
shows the structure of the hexagonal basal plane of UPd2Al3 
with the alignment of the a-axis (100) and 6-axis (-1,2,0), as 
used in the text. 



■^He cryostat set on a mechanical rotating stage at the 
top of a Dewar. By computer-controlling the two super- 
conducting magnets and rotating stage, we were able to 
rotate H with a misalignment of less than 0.5° from each 
axis, which we confirmed by the simultaneous measure- 
ments of the resistivity. 

In the present experiments, Kyy and Kzz were measured 
keeping H always perpendicular to q, i.e. the heat cur- 
rent perpendicular to the vortices dominates. This con- 
figuration (HJ-q) is important because when the angle 
between q and H is changed by rotating H, an addi- 
tional angular variation of the thermal conductivity ap- 
pears as a function of the angle between H and q due 
to the difference in the effective density of states (DOS) 
for quasiparticles traveling parallel to the vortices and 
for those moving in the perpendicular direction. In the 
field rotation experiments, the thermal conductivity was 
measured after field cooling above Tc- Consecutive mea- 
surements with an inverted direction produced little hys- 
teresis in the angular dependence of the thermal con- 
ductivity, which indicates that field trapping related to 
vortex pinning was negligible. 

Fi gure 1 depicts the temperature dependence of Kyy 
and Kzz in zero field. The magnitude and T-dependence 
of Kini and are similar to those reported in Ref. 
|24L liol. In UPd2Al3 there are three contributions that 
carry the heat; electrons, phonons and spin- waves. The 
spin- wave contribution appears to be negligible below T^, 
because a spin- wave has a finite gap of ~ 1.5 meV at the 
zone center and cannot contribute to the heat conduc- 
tion at low temperatures. The Wiedemann- Franz ratio 
L — at Tc is 0.95Lo for Kzz and is I.IGLq for Kyy, 
where Lq is the Lorentz number and p is the resistivity 
. These results indicate that the electron contribution 
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FIG. 2: Magnetic field dependence of the 6-axis thermal con- 
ductivity Kyy for H\\ a and H\\ c. Kyy is normalized by the 
normal state value just above the upper critical field Kyy. 



strongly dominates over the phonon contribution, at least 
below Tc- Unfortunately the present temperature range 
(the lowest temperature was 0.28 K) is not low enough to 
determine the detailed T-dependence and residual values 
of Kyy and Kzz, which would provide important informa- 
tion about the gap structure. We therefore do not discuss 
these subjects in this paper. 



Figure 2 shows the i7-dependence of Kyy for H 



yy 

and H\\ c below Tc- In both field directions, Kyy in- 
creases with H after an initial decrease at low fields. For 
H\\ c, Kyy increases almost linearly with H, Kyy oc H, 
at 0.36 K. The consequent minimum is much less pro- 
nounced at lower temperatures. As H approaches H^2 or 
H^2i i^yy shows a steep increase and attains its normal 
state value. In the normal state above Hc2, Kyy shows 
a slight decrease, which is scaled by the positive magne- 
toresisitance, supporting the conclusion that the electron 
contribution dominates the thermal conductivity. 

The i/-dependence of the thermal conductivity in the 
superconducting state is markedly different from that 
observed in ordinary s-wave superconductors, in which 
the thermal conductivity shows an exponential behav- 
ior with much slower growth with at <C i?c2 E3- 
The understanding of the heat transport for supercon- 
ductors with nodes h as p rogressed considerably during 
past few years ^3. S| • In contrast to classical super- 
conductors, the heat transport in nodal superconductors 
is dominated by contributions from delocalized quasipar- 
ticle states rather than bound states associated with vor- 
tex cores. The most remarkable effect on the thermal 
transport is the Doppler shift of the energy of quasipar- 
ticles with momentum p in the circulating supercurrent 
flow Vs [E(p) E{p) -Vs-p) miilltg. This effect 
(Volovik effect) becomes important at positions where 
the local energy gap becomes smaller than the Doppler 
shift term (A < Vg-p), which can be realized in the case of 
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superconductors with nodes. In the presence of hne nodes 
where the density of states (DOS) of electrons N{E) has a 
hnear energy dependence {N{E) cx E), N{H) increases 
in proportion to y/H. At high temperatures and low 
fields, where the condition y/H/H^ < T/Tc is satis- 
fied, the thermally excited quasiparticles dominate over 
the Dopplcr shifted quasiparticles. It has been shown 
that in this regime, while the Doppler shift enhances the 
DOS, it also leads to a suppression of both the impurity 
scattering time and Andreev scattering time off the vor- 
tices 32, 33^ ,43j ,44, 47}. This suppression can exceed 
the parallel rise in N{E) at high temperatures and low 
fields, which results in nonmonotonic field dependence 
of the thermal conductivity. As seen in other supercon- 
ductors with nodes, the initial decrease of the thermal 
conductivity disappears at low temperatures. The steep 
increase of iJc2, which is also observed in pure 

Nb, UPta, and Sr2Ru04, may be due to the enhancement 
of the quasiparticle mean free path caused by tunneling 
between vortex cores, which is possible near Hc2- Thus 
the if-dependence of Kyy in UPd2Al3, initial decrease 
at low field at high temperatures and linear behavior 
Kyy (X H at low temperatures, are in agreement with the 
existence of line nodes in A{H). We note that quite simi- 
lar ff-dependence has been reported in Sr2Ru04 with an 
anisotropic gap function |3^ . 

We here briefly discuss the band structure of UPd2 AI3. 
According to band calculations and de Haas-van Alphen 
measurements in the AF phase, the Fermi surface has 
four principle sheets of predominantly 5/-electron char- 
acter, which are labeled "cigar", " eggs" ," cylinder" , and 
"party hat" H Slill. The largest Fermi sheet with 
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FIG. 3: (a)Schematic figure of the gap structure with four 
line nodes perpendicular to the basal plane (vertical node). 
(b)Schematic diagram showing the regions on the Fermi sur- 
face that experience the Doppler shift in H within the basal 
plane. We have assumed dxy symmetry. <j!> — {H ,a) is the 
azimuthal angle measured from the a-axis. With H apphed 
along the antinodal directions, all four nodes contribute to 
the DOS, while for H applied parallel to the node directions, 
the Doppler shift vanishes at two of the nodes, (c) Four-fold 
oscillation of the DOS for H rotating in the basal plane. The 
DOS shows a maximum (minimum) when H is applied in the 
anitinodal (nodal) direction. 



large electron mass and the strongest 5f-admixture is the 
" cylinder" , which has the shape of a corrugated cylinder 
with hexagonal in-plane anisotropy. Therefore it is nat- 
ural to consider that the thermal conductivity is mainly 
determined by this Fermi sheet. In what follows, we as- 
sume that the cylindrical Fermi sheet with heavy mass is 
responsible for the superconductivity, neglecting all other 
sheets. 



III. ANGULAR-DEPENDENT 
MAGNETOTHERMAL CONDUCTIVITY 

Having established the predominant contribution of 
the delocalized quasiparticles in the thermal transport, 
the next issue is the nodal structure of UPd2Al3. An im- 
portant advantage of measuring the thermal conductivity 
is that it is a directional probe of the nodal structure. Re- 
cently measurement of the thermal conductivity or heat 
capacity with H applied various directions relative to the 
crystal axes has been established as a powerful method 
to determine the superconducting gap structures in k- 
space. At the heart of this method is the Volovik effect 
discussed previously. Since quasiparticles contribute to 
the DOS when the Doppler-shifted energy exceeds the 
local energy gap (A(fc) < Vg ■ p), an immediate conclu- 
sion is that the DOS depends sensitively on the angle 
between H and the direction of the nodes. The peri- 
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FIG. 4: (a) Schematic figure of the gap structure with line 
nodes parallel to the basal plane (horizontal node). Line 
nodes are located at the bottleneck (type-I) and at the AF 
zone boundary (type-II). Two line nodes are located at the 
bottleneck and the AF zone boundary (type-Ill) and at posi- 
tions shifted off the bottleneck (type-IV) . (b) Oscillations of 
the DOS for H rotating in the ac-plane for various gap func- 
tions. Two-fold oscillation with the same sign are expected 
for type-I, -II, and -III. On the other hand, for type-IV, an 
oscillation with a double minimum is expected. 
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odicity, phase, and shape of typical oscillations in the 
angular variation of the thermal conductivity give direct 
information of the gap structure, including the type of 
nodes (point or line) and their direction. Therefore, pos- 
sible allowed gap functions can be much more restricted 
in fc-spacc. It should be noted that this method does not 
rely on any specific mechanism of the superconductivity. 

We first discuss the case when H is rotated within 
the basal a6-plane. Here, for simplicity, we assume dxy 
symmetry with four line nodes located perpendicular 
to the basal plane, and assume a circular cross section 
of the Fermi sheet, neglecting the hexagonal in-plane 
anisotropy, as illustrated in Figs. 3 (a) and (b). Here 
(j)={H,a) is the azimuthal angle measured from the a- 
axis. One expects a significant anisotropy between the 
field induced DOS for H parallel to the nodal directions 
and H parallel to the anti-nodal direction. The Doppler 
effect should also lead to angular dependence of the DOS 
for H rotated within the basal plane. A minimum in 
the DOS occurs when H is along the nodal directions. 
Since, in that case, at two of the four nodes, the circu- 
lating currents are in the plane orthogonal to the mo- 
mentum vector at the nodes, the Doppler shift vanishes; 
only two nodes contribute to the DOS, as illustrated in 
Fig. 3. In contrast, when H is directed along the antin- 
odal directions, the Doppler shift is non-zero at some 
points in space for each node and all four nodes con- 
tribute to the DOS. As a result, the DOS oscillates with 
four-fold symmetry, as shown in Fig. 3(c). The ampli- 
tude of the DOS oscillation SN{E)/N{E) appears to be 
model dependent, but all calculations predict an ampli- 
tude ranging from 3% to 10% at zero temperature. In 
fact, a clear four- fold modulation of the thermal conduc- 
tivity, which reflects the angular position of the nodes, 
has been observed in the high- Tc cuprate YBa2Cu307_5 
[2^ I30I I31I I , heavy fermion CeCoIns [s^, , and organic 
K-(BEDT-TTF)2Cu(NCS)2 ^ with d-wave symmetries 
and in borocarbide YNi2B2C with s + g-wave symmety 
[3^ , while such an oscillation is absent in the fully 
gapped Y(Nio.95Pto. 05)2620 with s-wave symmetry [51| . 
demonstrating that thermal conductivity can be a rele- 
vant probe of the superconducting gap structure. In gen- 
eral, Kzz oscillates with n-fold symmetry corresponding 
to the number of vertical nodes n. 

We next discuss the situation when H is rotated within 
the ac-plane perpendicular to the basal a5-plane. The 
Fermi surface has open orbits along the c-axis. We here 
consider four gap functions in the AF Brillouin zone 
shown in Fig. 4(a). 

1. type-I : A horizontal node located at the bottleneck; 
A(A;) (X sin fc^c. 

2. type-II : A horizontal node located at the AF zone 
boundary; A(fc) cx cosfczC. 

3. type-Ill: A hybrid of type-I and -II. Two horizontal 
nodes located at the bottleneck and the AF zone 
boundary; A(A;) cx sin2kzC. 



4. type-IV : Two horizontal nodes located at posi- 
tions shifted off the bottleneck in the Brillouin zone; 
A(fc) cx cos 2kzC. 

The expected angular variations corresponding to these 
gap functions are shown schematically in Fig. 4(b). Here 
is the polar angle between H and the c-axis. The an- 
gular variation of the Doppler shifted DOS is a function 
of the relative angle between H and p. Therefore the 
twofold oscillations with the same phase are expected for 
type-I, -II , and -HI gap functions, in which the horizon- 
tal nodes are located at the position where p|| a6-plane; 
one cannot distinguish these three gap functions when 
the Fermi surface has an open orbit along the c-axis. For 
type-IV, one expects an oscillation with a double maxi- 
mum structure as a function of 9. 

We note that the angular variation is dominated by 
the Doppler shifted DOS at low fields. On the other 
hand, at high fields near Hc2, there are additional con- 
tributions which also cause the oscillation of the DOS, 
i.e. the anisotropy of the Fermi velocity and Hc2 j28ll5^ . 
Generally, the oscillation due to the Fermi velocity and 
due to Hc2 have the same sign. Therefore, it is desir- 
able to determine the nodal structure at low fields where 
the oscillation arising from Hc2 and the Fermi velocity is 
negligible. 



IV. RESULTS AND DISSCUSSION 

A. Vertical nodes 

We examine here the possibility of vertical line nodes 
perpendicular to the basal plane. Figure 5 shows 
K.zz{H^(j}) as a function of at 0.4 K, measured by ro- 
tating H within the basal plane. Above 0.5 T, a dis- 
tinct six- fold oscillation is observed in Kzz{H,(p), reflect- 
ing the hexagonal symmetry of the crystal. However, no 
discernible six- fold oscillation was observed below 0.5 T 
within our experimental resolution. Kzz can be decom- 
posed as Kzz = K^z + li^t where k^z ^ ^-independent 
term and n^f = Czf cos 6(j) has six- fold symmetry with 
respect to (f> -rotation. The sign of Cff is negative in 
the whole H range. The six-fold oscillation is clearly 
observed even in the normal state above Hc2- In con- 
trast, the amplitude of the sixfold oscillation, |Cf|'|/K"^, 
is less than 0.2% below 0.5 T, where k"^ is the c-axis 
thermal conductivity in the normal state just above Hc2- 
In Fig. 6, the ff-dependence of ICfl'l/K"^ is plotted as a 
function of H. In the region H^^ <H < H^^, \Czt\/f^zz 
is largest. This is because the sample is either in the 
normal or superconducting state with rotating H . With 
decreasing H, \Czf\/Kzz decreases and vanishes below 
0.5 T. 

We point out here that the AF magnetic domain struc- 
ture is responsible for the six-fold symmetry and the 
nodal structure is not related to the oscillation. For 
UPd2Al3, there are several AF ordered phases; three with 
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FIG. 6: Field dependence of the amplitude of the six-fold 
oscillation C^f normalized by the normal state value, k"^. 
H^2 = 3.2 T and ffc*2=3.12 T are the upper critical fields 
in H parallel to the a- and 6-axis, respectively. The inset 
shows the magnetic domain structures in parallel fields. The 
hexagon shows the basal plane. Populated magnetic domains 
with non-parallel spins are indicated by filled arrows. The 
open arrows indicate the antiparallel spins in the neighboring 
plane. The transition (i)-(ii) occurs at about 0.6 T. In the 
(i)-phase, the ordered moments point to the a-axis, forming 
domain structures with six-fold degeneracy. In the (ii)-phase, 
the domain structure changes to six-fold symmetry with re- 
spect to the iJ-rotation in the basal plane. 
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FIG. 5: Angular variation of the c-axis thermal conductiv- 
ity Kzz{H,(f)) normalized by the normal state value, k"^) a^t 
several fields at 0.4 K. H was rotated within the basal plane 
(see inset). A distinct six-fold oscillation is observed above 
0.5 T, while oscillation is absent at 0.5 and 0.3 T. At H=3.5 T 
(> Hc2), the system is in the normal state. 



the applied field in the easy (basal) plane, the (i)-,(ii)- 
, and (m)-phases, and one along the magnetically hard 
c-axis |5|. The magnetic structures in the (i)- and (ii)- 
phases are shown schematically in the inset of Fig. 6. The 
transitions, (i)-(ii) and (ii)-(iii), occur at 0.6 and 4.2 T, 
respectively, well below T/v . The fact that the sixfold os- 
cillation disappears in the vicinity the (i)-(ii) boundary 
strongly indicates that the six-fold oscillation is closely 
related to the AF magnetic structure. This idea is rein- 
forced by the fact that six-fold oscillation is observable 
even above Hc2- In the (i)-phase, the ordered moments 
point to the a-axis, forming domain structures with six- 
fold degeneracy inherent to the hexagonal crystal struc- 
ture. The spin structure is not effected by i?-rotation in 
the basal plane. On the other hand, for the (ii)-phase, 
if-rotation causes domain reorientation, giving rise to ro- 
tation of the ordered magnetic moments against a small 
anisotropy in the basal plane. As a result, the magnetic 
domain structure changes with six-fold symmetry with 



H rotation in the basal plane. Since the quasiparticles 
are scattered at the domain boundaries, k^z oscillates 
with six-fold symmetry in the (ii)-phase. Thus the six- 
fold symmetry observed in k^z above 0.5 T is most likely 
to be due to the magnetic domain structure. In addition, 
the in-plane anisotropy of Hc2 and the Fermi velocity is 
also important near Hc2- The most important finding 
is the absence of the oscillation in the (i)-phase where 
the spin structure is independent of H rotation. This 
definitely indicates that there are no nodes located per- 
pendicular to the basal plane, i.e. gap function in the 
basal plane is isotropic. 

B. Horizontal nodes 

We next examine the horizontal line nodes parallel to 
the basal plane. Figure 7 displays the angular variation 
of Kyy{H,9) for rotating H as a, function of 6 within 
the ac-plane at 0.4 K. A distinct oscillation with two- 
fold symmetry is observed in the superconducting state. 
In contrast to k^z, no discernible twofold oscillation was 
observed in the normal state above Hc2- i^yy{H,9) can 
be decomposed as Kyy = Kyy + Kyy where Kyy is a 6- 
independent term and Hyy = Cyy cos 26* is a term with 
twofold symmetry with respect to 6'-rotation. Figure 8(a) 
depicts the i/-dependence of Cyy at 0.4 K. For compar- 
ison, the i/-dependence of Kyy for H\\ c at T=0.36 K 
is plotted in Fig. 8(b). There are three regions denoted 
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FIG. 7: Angular variation of the &-axis thermal conductivity 
Kyy{H,6) normalized by the normal state value Hyy at several 
fields at 0.4 K. H was rotated within the ac-plane perpendic- 
ular to the basal plane (see inset). 
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FIG. 8: (a) //-dependence of the amplitude of the two-fold 
symmetry Cyy normalized by the normal state thermal con- 
ductivity Kyy at 0.4 K. The sign of Cyy is negative in the (I)- 
and (Ill)-regions, and is positive in the (Il)-region. (b) Field 
dependence of the fo-axis thermal conductivity Hyy for H\\ c 
at 0.36 K. 



(I), (II) and (HI), below Hc2- In the vicinity of Hc2 
((Ill)-region), where Kyy increases steeply with H, the 
sign of Cyy is negative and the amplitude \Cyy\/K,yy is of 
Here the normal state 



With decreasing H, Cyy changes sign 



the order of 10% 

just above Hc2- . , ^yy 

at about 2.3 T and become positive in the region where 
Kyy for H\\ c shows a linear _ff-dependence ((Il)-region). 
Below about 0.25 T, where the second sign change takes 
place, Kyy decreases with H ((I)-region). In this region, 
Cyy becomes negative. 

We here address the origins of the observed two-fold 
oscillation. The disappearance of the oscillation above 
Hc2, together with the fact that there is only one mag- 
netic phase in this configuration completely rule out 
the possibility that the origin is due to the magnetic do- 
main structure. There are two possible origins for the 
oscillation; the nodal structure and the anisotropy of the 
Fermi velocity and Hc2- Obviously, as discussed pre- 
viously, a large two-fold oscillation with negative sign 
observed in the (Ill)-region arises from the anisotropies 
of the Fermi velocity and Hc2- This immediately indi- 
cates that the two-fold symmetry with positive sign in 
the (Il)-region originates not from these anisotropies hut 
from the quasiparticle structure associated with the nodal 

-i20 



the (Il)-region is a few percent, which is quantitatively 
consistent with the prediction considering the Doppler 
shifted DOS. We also note that the second sign change 
at low fields in the (I)-region is compatible with the nodal 
structure. In this region, as discussed previously, the H- 
dependence of the thermal conductivity is governed by 
the suppression of th e quasip article scattering rate. As 
discussed in Ref. [H |H iH ED, El El 113 > the oscilla- 
tion arising from the anisotropic carrier scattering time 
associated with the nodal structure also gives rise to the 
oscillation of Kyy{H,d) as a function of 0. In this case the 
sign of the oscillation is opposite to that arising form the 
Doppler shifted DOS in the (Il)-region. These consider- 
ations lead us to conclude that UPd2Al3 has horizontal 
nodes. In addition, the fact that there is a single maxi- 
mum structure in the angular variation of Kyy {H,0) indi- 
cates that horizontal line nodes are located at positions 
where the condition p\\ ab in the Brillouin zone is satis- 
fied. Thus the allowed positions of the horizontal nodes 
are restricted at the bottleneck and AF zone boundary. 



V. CONCLUSION 

For comparison, the angular variations of Kyy and Kzz 
gap function. In addition, the amplitude of Cyy/ Kyy in at low fields are shown in Fig. 9. While the amplitude 
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FIG. 9: Angular variation of the 6-axis thermal conductiv- 
ity K,yy/i-iyy with rotating H within the ac-plane and of the 
c-axis thermal conductivity K^t/i^zz with rotating H within 
the basal afe-plane at r=0.4 K and at H=Q.5 T. The ampli- 
tude of the six-fold oscillation in K®f / kJ^ is less than 0.2% if it 
exists. Inset: Schematic figure of the gap function of UPd2 AI3 
determined by angle resolved magnetothermal transport mea- 
surements. The thick solid lines indicate horizontal nodes 
located at the AF zone boundaries. 



of the two-fold oscillation Cyy/Kyy is 3%, which is quan- 
titatively consistent with the Doppler shifted DOS, the 
amplitude of the six-fold oscillation f / Kzz is less than 
0.2%, which is more than 20 times smaller than the os- 
cillation expected from the Doppler shifted DOS in the 
presence of nodes. Combining the results, we arrive at 
the conclusion that the gap function is isotropic in the 
basal plane and has horizontal node. The order parame- 
ters allowed, by thermal conductivity measurements, are, 

1. A(fc) = Aosin kzC, 

2. A(fc) = Aosin2A;2C and 

3. A(fc) = AocosfczC. 

which are shown in the type-I, -II and -III gap struc- 
tures in Fig. 4(a). The first and second represent spin 
triplet gap functions, and the third represents a spin 
singlet gap function. As discussed previously, the i?c2, 
NMR, and /.tSR experiments all favor spin singlet pair- 
ing. We note that spin singlet pairing is also supported by 
tunneling spectroscopy measurements. In anisotropic su- 
perconductors, zero-bias conductance peak appears due 
to mid gap Andreev states where the injected and re- 
flected electrons feel the different signs of the pair po- 
tentials [53. In the experiment reported in Ref. [lol |. 
in which tunneling conductance is measured along the 



c-axis, a zero-bias conductance peak appears for triplet 
paring while it is absent for singlet pairing. The ab- 
sence of the zero-bias conductance peak supports singlet 
pairing. Furthermore, the third gap function is compati- 
ble with the constraint implied by the neutron resonance 
peak, A(fc) = — A(fc + Qo)- These considerations lead 
us to conclude that the gap function of UPd2Aly, is most 
likely to be A(fe) = Aocosfc^c with d-wave symmetry, 
shown in the inset of Fig. 9. 

Theoretical models with various nodal structures for 
UPd2Al3 have been proposed. The determined gap 
structure cannot be explained in the framework of the 
electron-phonon coupling mechanism. There are two 
classes of theories of the mechanism of the unconventional 
superconductivity in UPd2Al3. The first is superconduc- 
tivity mediated by the exchange of AF spin- wave excita- 
tion (magnetic excitons); the coupling between localized- 
and itinerant-subsystem is essential for the occurrence 
of superconductivity 0, 0, I3- The second is a mech- 
anism based on the AF spin fluctuation arising from 
the strong electron-electron correlation in the itinerant 
subsystem, ignoring the localized subsystem |54j . This 
model predicts vertical nodes, which is inconsistent with 
the present results. The existence of horizontal nodes 
strongly supports the interpretation inferred from the 
resonance peak observed at the AF Bragg point in the 
neutron inelastic scattering experiments. The horizontal 
node located at the AF zone center indicates that pair 
partners cannot reside in the same basal plane. The in- 
terlayer pairing appears to indicate that strong dispersion 
of the magnetic excitation along kz causes the pairing, as 
suggested in the magnetic exciton mediated superconduc- 
tivity model 0, 0j I2] ■ This dispersion is connected with 
the fluctuation of the out-of-plane component of the or- 
dered moments, because the in-plane moments are satu- 
rated at Tc ^ T/v and should not have large fluctuations. 
The isotropic gap function in the basal plane implies that 
the pairing interaction in the neighboring planes strongly 
dominates over the interaction in the same plane. This 
is incompatible with the theory in which only the elec- 
tron correlation in the basal plane is responsible for the 
superconductivity. Although the pairing interaction in- 
ferred from the determined gap function should be fur- 
ther scrutinized, the recent results imply that the inter- 
layer pairing interaction associated with the AF interac- 
tion is most likely to be the origin of the unconventional 
superconductivity in UPd2Al3. 
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